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Bilayer graphene is a highly tunable material: not only can one tune the Fermi energy using 
standard gates, as in single-layer graphene, but the band structure can also be modified by external 
perturbations such as transverse electric fields or strain. We review the theoretical basics of the 
band structure of bilayer graphene and study the evolution of the band structure under the influence 
of these two external parameters. We highlight their key role concerning the ease to experimentally 
probe the presence of a Lifshitz transition, which consists in a change of Fermi contour topology as a 
function of energy close to the edges of the conduction and valence bands. Using a device geometry 
that allows the application of exceptionally high displacement fields, we then illustrate in detail the 
way to probe the topology changes experimentally using quantum Hall effect measurements in a 
gapped bilayer graphene system. 


INTRODUCTION 

The shape of the Fermi surface is crucial for under¬ 
standing the electronic properties of metals [1]. As first 
noticed by Lifshitz [2] , changes in the Fermi surface topol¬ 
ogy cause anomalous behavior of thermodynamic, trans¬ 
port and elastic properties of materials [3]. Intuitively, 
the simplest way to observe such an electronic topological 
transition, also known as a Lifshitz transition, is by tun¬ 
ing the Fermi level to the singular point in the band struc¬ 
ture where the change of topology takes place. However, 
this usually requires considerable variations of the elec¬ 
tron density that for bulk materials can only be achieved 
by the experimentally inconvenient preparation of nu¬ 
merous samples with different chemical compositions [4- 
9]. Alternatively, the Fermi surface can be deformed by 
exposing the sample to extreme conditions like high pres¬ 
sures [2, 10, 11] or strong magnetic fields [12, 13]. 

In this article we discuss the theoretical prediction and 
an experimental study of the Lifshitz transition in bi¬ 
layer graphene [14], two coupled layers of carbon atoms 
arranged in a hexagonal lattice [Fig. 1(a)]. Its elec¬ 
tronic dispersion relation contains saddle points, simi¬ 
lar to those well known from graphite [15-18]. However, 
as opposed to bulk graphite, external perturbations like 
strain [19, 20] or electric fields [21] can modify the topol¬ 
ogy of the electronic dispersion and change the energy of 
the Lifshitz transition which connects regions of differ¬ 
ent Fermi contour topologies [Fig. l(c-e)]. Moreover, due 
to the two-dimensional nature of bilayer graphene, its 
chemical potential and its topology can be tuned with 
electrostatic gates [22], simplifying experimental studies 
of the Lifshitz transition. 

We first discuss in detail the band structure of bi- 
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FIG. 1. (a) Schematic lattice structure of bilayer graphene. 
Also shown are all the electron hoppings mentioned in this 
work, (b) Electronic structure of bilayer graphene on the en¬ 
ergy scale of 3 eV from the neutrality point (marked by the po¬ 
sition of the black hexagon which shows bilayer graphene Bril- 
louin zone). The blue (yellow) surfaces show the low-energy 
(split) bands. Red circles indicate the K and K' valleys, (c-e) 
Low-energy electronic structure at the K valley of an (c) un¬ 
perturbed, (d) homogeneously strained and e) gapped bilayer 
graphene. Grey discs cut through the dispersion at the energy 
of the Lifshitz transition and the corresponding cross-section 
is shown in the right of each of the panels. 


layer graphene in the absence and presence of an exter¬ 
nal quantizing magnetic field perpendicular to the car¬ 
bon layers. We follow on with the analysis of the ef¬ 
fects on the electronic spectrum of external parameters, 
such as homogeneous strain and electric fields applied 
across the layers. We then describe state-of-the-art bi¬ 
layer graphene devices and the experimental signatures 
of the Lifshitz transition in this material. 
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I. ELECTRONIC BAND STRUCTURE OF 
NEUTRAL BILAYER GRAPHENE 


A. In the absence of the magnetic field 


Bilayer graphene, shown in Fig. 1(a), is composed of 
two graphene [23] layers, that is, two layers of hexag- 
onally arrayed carbon atoms. The unit cell consists of 
four atoms which we label as A 1 , B 1 , A2 and B 2, with 
A and B labeling the sublattice within the same layer 
and the index 1 (2) denoting the bottom (top) layer. 
In the present study, we focus on Bernal-stacked bi¬ 
layer graphene, i.e. atom A2 lies on top of atom Bl. 
For the electronic properties, it is enough to consider 
only the delocalized 7r-electrons, one per carbon atom. 
Within the tight-binding approach, the essential features 
of the electronic dispersion are captured by taking into 
account the in-plane nearest neighbor electron hopping 
7o = Iaibi = 1 A 2 B 2 ~ 3 eV, the vertical interlayer 
coupling 71 = r )A 2 Bi ~ 0.4 eV and the skew inter¬ 
layer hopping 73 = r )AiB 2 ~ 0.4 eV [24]. The result¬ 
ing electronic bands are shown in Fig. 1(b) within the 
hexagonal Brillouin zone containing two inequivalent cor¬ 
ners K and A', also called valleys. Two of the bands, 
which are shown in blue in Fig. 1(b) and which we re¬ 
fer to as low-energy ones, touch at the center of the 
valley at the energy which marks the neutrality point 
- the position of the chemical potential in neutral bi¬ 
layer graphene. The other two, called split bands and 
shown in yellow in Fig. 1(b), are split off the neutrality 
point by the energy e ~ ±71. In the vicinity of the val¬ 
ley, the electronic band structure can be described by a 
four-band Hamiltonian [14] written in the basis of sublat¬ 
tice Bloch states ((j)Ai,4>B2,4>A2,4>Bi) T in valley K and 
(0B2, 4>ai,4>biAa2) T in valley K ', 

0 V 37 f 0 

f T Vs7T^ 0 V7 f 

rl — , 

0 VTZ' 0 
_ V7T 0 £71 

Above, £ = 1 (£ = — 1) refers to valley K ( K' ), v = 

Q ^| 7 ° , v 3 = Q ^| 73 and 7 r = p x -\- ip y with the momentum 
p = ( Px,Py ) measured from the center of the valley. 

Using Hamiltonian H in Eq. (1), the dispersion of the 
low-energy bands can be written as 


f7T' 
0 

£71 

0 


(1) 


e “ = “71 bl + 2«V + v Ip 2 - [li + 4 7i v 2 P 2 (2) 

~ v 3P 2 (^7l - 4u V - v 3P 2 ) 

1 

1 

+ 8 £ 7 i^ 3 v 2 p 3 cos3(^] 2 

where a = 1 (<a = — 1 ) denotes the conduction (va¬ 
lence) band and p = \p\. While for a single-layer of 
graphene, the in-plane coupling 70 leads to the linear 



Dirac-like dispersion of electrons [25] with the “effective 
speed of light” v « 10 6 m/s [23], the coupling 71 between 
the two layers turns the linear dispersion into a quasi- 
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parabolic one, giving e ~ ±^, at intermediate energies 
(71 (^) 2 < H < 71 ) [14]. The effective mass is given by 
m = 7^/2 ~ 0.034m e , where m e is the free electron mass. 
The skew interlayer hopping 73 introduces trigonal warp¬ 
ing of the electronic dispersion, visible in the shape of the 
constant-energy contour in Fig. 1(c). Most importantly, 
in unperturbed bilayer graphene it leads at low energies 
(|e| < ^ 7 i(^f ) 2 = enr ~ 1 meV) to the fragmentation of 
the singly-connected Fermi line into four separate pieces 
[14], as shown in the right of Fig. 1(c). This is because 
the band structure at the neutrality point consists of four 
(anisotropic) mini Dirac cones. Because the energy of the 
Lifshitz transition 6 lt at which this fragmentation of the 
Fermi line occurs is very close to the neutrality point, 
the carrier density required to shift the chemical poten¬ 
tial to clt is small, tilt ~ 2 x 10 10 cm -2 . In fact, this 
carrier density is smaller than or comparable to the local 
charge density fluctuations induced by charged impuri¬ 
ties in most of the substrates used for graphene. For that 
reason, as described in Sec. II, extremely clean flakes are 
required to probe the Lifshitz transition in unperturbed 
bilayer graphene. 

The trigonal warping that leads to the Lifshitz transi¬ 
tion in bilayer graphene results in a similar spectral fea¬ 
ture in the vicinity of the A-point of the Brillouin zone 
of bulk graphite [15, 18]. Because the chemical poten¬ 
tial in this case, like for unperturbed bilayer graphene, is 
close to the Lifshitz transition, signatures of the change 
of topology of the band structure can be observed for ex¬ 
ample in cyclotron resonance experiments [18]. However, 
despite the similarities between the band structures of 
graphite and bilayer graphene, we emphasize that the for¬ 
mer lacks the tunability of the latter, both with regards 
to changing the Fermi level and modifying the Lifshitz 
transition with external factors. 


B. Electronic states in the presence of a strong 
perpendicular magnetic field 

In a sufficiently strong magnetic field perpendicular to 
the graphene layers, the electronic band structure de¬ 
scribed before undergoes Landau quantization and forms 
a ladder of discrete energy levels. We describe the theo¬ 
retical procedure used to calculate the Landau level (LL) 
structures shown in this paper in Appendix A. In Fig. 
2 we show a comparison between the LL spectra calcu¬ 
lated using the Hamiltonian in Eq. (1), for 73 =0 (solid 
red lines) and 73 7 ^ 0 (solid black lines). At high mag¬ 
netic fields the two spectra are essentially identical, with 
the energy of a LL n a (a is used to distinguish between 
the electron and hole levels with the same n) described 
by e n a = ah^^/n(n — 1) [14]. According to this for¬ 
mula, Landau levels n = 0 and n = 1 have the same 
energy, e = 0, so that the zero-energy LL, characteristic 
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FIG. 2 . Comparison of the Landau level spectra calculated 
for gapless bilayer graphene without taking into account the 
73 coupling (solid red lines) and with 73 included (solid black 
lines). For the latter case, we assumed vs/v = 0.1. 

for graphene, exhibits an extra double orbital degeneracy. 
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This additional degeneracy results in a step of 8 ^ in the 
measured Hall conductance when the Fermi level crosses 
the neutrality point and suggests that the filling factors 
v = ±4 display the largest activation gap [14, 26]. The 
higher LLs form an almost equidistant staircase of elec¬ 
tronic states and are fourfold degenerate due to the valley 
and spin degeneracies. In high quality bilayer graphene, 
strong electron-electron interaction can however lift these 
degeneracies, as will be discussed later. 

As the low-energy band structure contains four mini 
Dirac cones, at very low magnetic fields the zero-energy 
Landau level has an additional 16-fold degeneracy (four 
Dirac cones x2 valleys x 2 spins). However, already at 
magnetic fields 5 ^ 0.1 T, the inverse of the magnetic 

length Xb = \J~^ is comparable with the distance in mo¬ 
mentum space between the mini Dirac points. This leads 
to magnetic breakdown [27, 28] as the separate pieces of 
the Fermi line are no longer resolved individually. As a 
result, the 16-fold degeneracy of the zero-energy Landau 
level is lifted at B « 0.1 T and the high-field structure 
described above is formed. Due to the presence of the 
Lifshitz transition, in high quality samples one might ex¬ 
pect filling factors v = ±8 to be the most robust at the 
lowest magnetic fields. 

In this theoretical part, we described the role of each 
interatomic coupling constant on the band structure and 
revealed the importance of the skew interlayer coupling 
73 , which gives rise to a Lifshitz transition close to the 
charge neutrality point. We showed that observing this 
transition is challenging because it takes place at low en¬ 
ergies compared to the usual amount of disorder present 
in bilayer graphene samples. Obtaining high quality sam¬ 
ples is therefore one key requirement to experimentally 


access this phenomenon. In the next section, we will ex¬ 
plain how high quality graphene can be obtained in ex¬ 
periments and how its band structure can be influenced 
by external parameters. 


II. THE TUN ABILITY OF BILAYER 
GRAPHENE 

Since its discovery, graphene has raised a lot of experi¬ 
mental interest, both from the fundamental point of view 
(offering, among other things, access to relativistic phe¬ 
nomena in a solid state system) and because of potential 
electronic applications [23]. However, graphene quickly 
appeared to be only that “perfect” in the theory world: 
graphene, once isolated from graphite flakes, has to be 
deposited on a substrate and further processed to be ex¬ 
perimentally probed. These fabrication steps necessarily 
introduce disorder. In this part, we focus on explaining 
how the quality of graphene devices can be improved by 
finding substrates alternatives to Si/SiO 2 and how this 
allows to experimentally access the Lifshitz transition. 


A. Substrate-induced disorder in single and bilayer 
graphene 

In early graphene-based transport experiments, 
graphene was deposited on Si/Si02 substrates [22, 29, 
30], which offer both good optical contrast to locate 
graphene and a way to electrostatically tune the Fermi 
energy using doped silicon as a backgate. However, it was 
rapidly found that this substrate, as convenient as it is, is 
a non-negligible source of disorder: the surface is rough, 
contains impurities, charge traps as well as charged sur¬ 
face states, which all limit the electronic performance of 
the resulting devices [31-34]. To circumvent this issue, 
two main paths have been investigated. 

The first one consists in suspending graphene to re¬ 
move the substrate altogether. Until recently, this 
yielded the best quality graphene ever achieved [35, 36]. 
This technique, however, results in fragile devices with 
complicated architectures if additional electric gates are 
required. It also limits the range of vertical electric fields 
that one can apply to graphene before it collapses. At 
the same time, this approach could potentially be used to 
induce controlled strain in the suspended graphene mem¬ 
brane [37]. As discussed in the next Section, this leads to 
interesting consequences for the band structure of bilayer 
graphene. 

In parallel, another option was developed by Dean et 
al. [38]: using hexagonal boron nitride (h-BN) as a sub¬ 
strate. This material is both a good substrate (low rough¬ 
ness and free of charge traps) and a good dielectric (it 
is very robust and allows the application of high elec¬ 
tric fields). The first process developed by Dean et al., 
called the “dry transfer technique”, showed a clear en¬ 
hancement of the charge carrier mobilities [38]. However, 
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these numbers could still not compete with suspended 
devices, mostly because of the difficulty to get rid of 
process-induced residues. But very recently, the same 
group developed a new technique [39], simultaneously 
with other groups [40], called the “pick-up” technique, 
allowing for one-dimensional contacts. This new process 
allows one to create h-BN/Graphene/h-BN stacks that 
can be contacted at the edge, leaving the graphene free 
of any polymer residue. This also allows the graphene-on- 
boron nitride technology to now compete with the quality 
of suspended devices, with the advantage of being free of 
strain and to allow the application of large gate voltages 
through the flake. 

Another advantage of the encapsulation of bilayer 
graphene lies in the fact that one can easily pattern top 
gates on the surface of the top h-BN: together with the 
Si-backgate, they allow for the application of a poten¬ 
tial asymmetry between the top and bottom layer of the 
bilayer flake. Such an asymmetry has profound conse¬ 
quences for the band structure of bilayer graphene and, 
as discussed later in this review, on the Lifshitz transition 
present in this material. 

We should however highlight at this stage that the 
substrate removal or replacement does not solve all the 
problems identified by transport experiments: besides 
the clear improvement of the bulk transport quality, the 
problems concerning nanostructures remain unchanged 
[41], as it seems that the performance of these devices is 
mostly dominated by the edges [42]. 

In the following, we will explain in more detail the 
theoretically expected effects on the band structure of 
strain and of the field-induced asymmetry, and we will 
show their influence on the Lifshitz transition. 


B. Influence of strain 

In this section, in particular the influence of (homo¬ 
geneous) strain on the topology of the low-energy bands 
and the Lifshitz transition energy clt- Strain can poten¬ 
tially be induced in suspended flakes, either intentionally 
[37] or unintentionally, for example as a side effect of 
the fabrication process. We neglect here the out-of-plane 
deformations of the graphene sheets and consider an un¬ 
deformed bilayer graphene crystal, shown in Fig. 3(a), 
which undergoes homogeneous deformation characterized 
by an angle 6 between the principal axis of strain tensor 
and the x-axis of the coordinate system as depicted in 
Fig. 3(b), and its eigenvalues, 8 and 8 '. Note that after 
such a deformation, the atomic sites A2 and Bl remain 
on top of each other. We further include shear deforma¬ 
tions by allowing a shift of the top layer with respect to 
the bottom one, by dr = (Sx,5y), as shown in Fig. 3(c). 
As a result of homogeneous strain, (1) the Brillouin zone 
is distorted due to the real space deformation of the lat¬ 
tice, (2) distances between carbon atoms change imply¬ 
ing modification of the hopping integrals. The former 
is a purely geometrical effect leading to a shift of the 
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FIG. 3. Top view of (a) unstrained, (b) strained, (c) strained 
and displaced graphene layers in bilayer graphene. 


valleys from their position in the reciprocal space. The 
latter means that the cancellation between the contribu¬ 
tions of the three 70 ( 73 ) hops for the coupling of plane 
wave Bloch states on A1/A2 and B1/B2 sublattices in 
the (shifted) valleys K and K' of the Brillouin zone is no 
more exact, leading to the Hamiltonian [19, 20] 

0 V^TT + W 0 V7T^ 

+ re* 0 vtt 0 /g\ 

0 V7V^ 0 £71 

V7T 0 £71 0 

where 

W = ^(r /3 - r]o)^ 3 e~ %20 (5 - S') - ^73^3^——, tany> 

Above, tab is the equilibrium carbon-carbon distance 
and % = and % = quantify the 

change of the intralayer A-B and interlayer A1-B2 hop¬ 
pings, respectively, upon the change of the distance be¬ 
tween the coupled carbon atoms. Note that w is in gen¬ 
eral a complex number with the phase specified by the 
directions of the principal axis of strain and the shear 
deformation. In the case of the former, rotation by i r 
leaves w unchanged (it does not matter by which of the 
opposite edges the material is stretched). Also, for the 
same eigenvalues of the strain tensor and ignoring shear, 
applying strain along the zigzag or armchair directions 
is equivalent to switching the sign of w. For shear de¬ 
formation, note that to simplify notation the angle ip is 
defined as the one between the negative direction of the 
y axis and the direction of shear. 

Based on the analysis of Raman spectra on monolayer 
graphene, it has been estimated [19] that 1% of strain 
can lead to w ~ 6 meV. In the same work it has also 
been established that electron-electron interaction would 
further enhance the strain-induced effects. This can be 
compared to the Lifshitz transition energy 6 lt ~ 1 meV 
in the unperturbed bilayer graphene. 

The resulting low-energy electronic dispersions of 
strained bilayer graphene in the vicinity of the valley K 
are displayed for representative values of (JNn, Ssw) in Fig. 
4. The dispersion for the case of the unperturbed bilayer 
graphene, that is for w = 0 , is shown in the center of 

the top row in Fig. 4. As mentioned before, it is quasi- 

2 

parabolic at energies e clt = —^ > with v 3 responsible 
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for the trigonal warping of isoenergetic lines. At e = 6 lt, 
the isoenergetic line undergoes a Lifshitz transition and 
splits into four Dirac cones [14]. The remaining disper¬ 
sions in the top row depict the band structure in the 
vicinity of the neutrality point for $lw ^ 0 , ^sw = 0 . As 
the value of $lw is increased (right side of the top row 
in Fig. 4), the two side cones positioned off the p x axis 
and the central cone move closer, as shown in the graph 
for w = 6 lt- Those three cones collide for w = 36lt, 
and for w > 3clt only two Dirac cones remain, as shown 
in Fig. 4 for w = 5clt- If Jim < 0 is negative (left side 
of the top row in Fig. 4), the central cone and the side 
cone positioned on the x axis approach each other and 
collide for w = — 6 lt, creating a local, quasi-parabolic 
minimum. This minimum persists for some range of the 
strain, although it lifts off the e = 0 plane, as shown in 
the graph for w = — 5clt- Eventually, for $lw = — 96lt, 
the minimum merges with a saddle point and with further 
decrease of again only two Dirac cones remain in the 
spectrum, as shown in the spectrum for w = —9.5clt- A 
contrasting situation of = 0 , ^ 0 , is presented in 

the graphs in the column in the middle of Fig. 4 (due to 
symmetry, dispersions for opposite values of are mir¬ 
ror reflections of each other with respect to the p y = 0 
plane and hence only situation of ^sw < 0 is described 
here). Again, starting from the unperturbed system and 
increasing the magnitude of leads to two of the cones 
moving closer to each other and colliding with a creation 
of a local, quasi-parabolic minimum which lifts off the 
e = 0 plane (see the dispersion for w = — i3eirr) and even¬ 
tually disappears, leaving only two Dirac cones (graph 
for w = —z56lt)- Three qualitatively different regimes in 
the (JRw, Ssw) space can be identified: ( 1 ) the spectrum 
contains four Dirac cones, ( 2 ) the spectrum contains two 
Dirac cones and a local minimum, created by a collision 
of two of the cones and which in general does not touch 
the 6 = 0 plane, and (3) the spectrum contains only two 
Dirac cones. The extent of those regimes is shown in 
the bottom left of Fig. 4, together with the points in 
this space corresponding to spectra discussed above and 
shown in Fig. 4. Note that symmetry of the lattice re¬ 
quires that, for no shear, rotating the principal axes of 
strain by 7 r /3 leads to the same deformation. In the case 
of pure interlayer shear, choosing directions related by 
27r/3 rotation also results in the same physical situation. 
Both are reflected by the 2tt/3 symmetry of the graph in 
the bottom left in Fig. 4 (as the angle 6 of applied strain 
appears doubled in the definition of w). 

The strain-induced deformation of the low-energy elec¬ 
tronic dispersion gives rise to distinctive features in the 
low-energy density of states (DoS) in each of the strain 
regimes. In particular, existence of a local parabolic 
minimum results in a sharp step in the DoS, while sad¬ 
dle points lead to van Hove singularities. Note that for 
w 7 ^ 0 , the symmetry between the side cones is bro¬ 
ken and the saddle points are no longer all found at the 
same energy 6 lt> leading to more than one van Hove 
singularity and more than one Lifshitz transition. Densi¬ 


ties of states representative of all the regimes are shown 
in the bottom right of Fig. 4 for re = 0, w = —i3eLT 
and w = — i5eLT, that is those of the points marked in 
red in the diagram in the bottom left of Fig. 4 that he 
on the 3?re = 0 line. Indeed, for re = 0 only one van 
Hove singularity is present, whereas two can be seen for 
re = —i3eLT- A sharp step, corresponding to the contri¬ 
bution of a quasi-parabolic part of the band structure to 
the density of states, is also visible for the latter case. 
Finally, for re = — i5eur, the sharp step is washed out (as 
the local minimum merged with the saddle point) and 
only one van Hove singularity remains. 

The transformation of electron dispersion by homo¬ 
geneous strain leads to the modification of the bilayer 
graphene LL spectrum. The examples of numerically 
calculated LLs are shown for the electronic dispersion 
in the top row of Fig. 4 for low magnetic fields, B < IT. 
Both for small and large strain, the high-magnetic-field 
end of the LL fan plot can be described neglecting the 
trigonal warping, as mentioned in Section IB. At low 
fields, such that huj c (B) < rar>f, the largest gap in the 
LL spectrum is between the e = 0 and next excited LL, 
suggesting the persistence of filling factor v = ±8 in the 
quantum Hall effect (QHE). After strain causes the an¬ 
nihilation of two out of four Dirac points, the e = 0 level 
becomes eightfold degenerate, and, for strain \w\ 6 lt, 
only filling factors v = +4 and u = — 4 persist in the 
low-field QHE in bilayer graphene: the largest energy 
gap in the LL spectra is between the eightfold degener¬ 
ate level at e = 0 and next excited level, whereas the 
rest of the spectrum is quite dense. This eightfold degen¬ 
eracy is topologically protected and it also appears in a 
rotationally twisted two-layer stack [43] . Importantly, for 
intermediate strains such that the low-energy electronic 
dispersion displays a quasi-parabolic local minimum, the 
LL structure contains a Landau level origins of which 
can be traced to electronic states in that minimum - its 
dispersion as a function of the magnetic field B is linear 
and its energy with vanishing B does not go to zero but 
rather follows the position of said minimum. 


C. Influence of the interlayer asymmetry u 

As mentioned earlier, the encapsulation of bilayer 
graphene enables the use of top gates. This allows ap¬ 
plying a high field-induced asymmetry between the two 
layers which strongly influences the band structure.. 


1. Inducing an asymmetry u with 73 = 0 : the Mexican hat 

One of the most interesting things about bilayer 
graphene as a material is the fact that, since it possesses 
two layers, one can try to address each layer indepen¬ 
dently. To do so, one can tune the bottom layer to a po¬ 
tential u/2 and the top layer to a potential —u/ 2. This 
way, an asymmetry of amplitude u is induced between 
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FIG. 4. The low-energy band structures for strained bilayer graphene. For the band structures shown in the horizontal row, 
the corresponding Landau level fans are shown underneath. Densities of states calculated for the electronic spectra displayed in 
the vertical column in the center are shown in the bottom right. Each electronic spectrum is assigned a label colored according 
to the characteristic topology of the spectrum. The extent of the parametric regimes of complex w, see Eq. (3), distinguishing 
between three characteristic topologies of the strained bilayer graphene spectrum is displayed using those colors in the bottom 
left. In this graph, red dots show the points in the parametric space corresponding to the spectra presented in this figure. 


the two layers. The resulting Hamiltonian is now given 
by 


2 U 

0 

0 

V7T^ 


0 

_£ 

V7T 

0 

(4) 

0 

4 

V7T 

z 

2 

£71 

VIT 

0 

f7i 

! u - 



To induce this asymmetry, chemical doping can be used 
[44] or external gates [45] can be patterned. The relation 
between u and the applied voltages has been investigated 
in Ref. [46, 47]. This induced asymmetry opens a band 
gap, as shown in Fig. 5(a) [14, 46, 48, 49]. The gap gives 
rise to a characteristic shape of the conduction/valence 
band edge, a so called “Mexican hat”. This is shown 
in Fig. 5(b). For an induced asymmetry u = 100 meV, 
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FIG. 5. (a) Effect of the interlayer asymmetry on the low 
energy band structure of bilayer graphene in the proximity of 
the A-point, calculated taking into account only the in-plane 
coupling 70 and the interlayer coupling between dimers 71 
(73 = 0): a band gap is opened (here u is set to 100 meV). (b) 
Zoom at the top of the valence band: the so-called “Mexican 
hat” is visible. For u = 100 meV, the amplitude of the hat is 
A u — 1.6 meV. 

the amplitude of the central dip is A u = 1.6 meV. The 
corresponding density of states is shown in Fig. 6 (black 
dashed curve). 

This asymmetry also has consequences for the mag¬ 
netic field behavior, as it lifts the valley degeneracy of the 
Landau levels [14] . This is because the Landau level wave 
function is not distributed equally among the layers, with 
the most striking example being that of the zero-energy 
Landau level. For this level, electronic states in the valley 
K are formed by orbitals that sit mostly (and exclusively 
in the case of gapless bilayer) on the Al-sites, while in the 
valley A', these states occupy dominantly the B 2 sublat¬ 
tice. Thus, since for nonzero interlayer asymmetry these 
sites are set to different potentials, electronic states in A 
are no longer degenerate with states in K '. Such lifting of 
the valley degeneracy occurs for all Landau levels but the 
magnitude of the energy splitting decreases with increas¬ 
ing Landau level index [14]. In experiments, even though 
the lifting of all the degeneracies has been observed for 
the lowest LL, temperature and energy broadening usu¬ 
ally hinder the splitting for higher LLs. 



FIG. 6 . Density of state calculated for the case with trigo¬ 
nal warping included (red curve) and without (73 = 0 , black 
dashed curve). The height of the peaks in the DoS was lim¬ 
ited for visual purposes to enable better comparison between 
other features in the graph. 


2. Influence of 73 in a gapped bilayer graphene system 

As mentioned previously, accessing experimentally the 
Lifshitz transition in bilayer graphene is challenging, 
ffowever, by using the possibility of applying a vertical 
electric field, one can realize such a situation experimen¬ 
tally. 

Considering the following electron-hole symmetric 
Hamiltonian, 


’ 2 U 

^ 3 7f 

0 

V7V^ 

~ t 

V 3 7T ] 

£ 

2 

V7T 

0 

0 

V7V' 

£ 

2 

C7i 

_ VI f 

0 

£71 

§M_ 


one can now set the interlayer asymmetry to u = 
100 meV and observe the consequences on the band struc¬ 
ture. As shown in Fig. 7(a), the band gap is open in a 
similar way as in the case where no skew interlayer hop¬ 
ping was included. However, already on this scale of 
±0.8 meV, one can see that the trigonal perturbation is 
still present. Zooming at the top of the valence band, as 
done in Fig. 7(b), we see that the band gap opening has a 
strong effect on the previous four Dirac cones: only three 
peaks are left and they do not exhibit a linear dispersion 
anymore. The central peak has been pushed down by the 
71 parameter, by the same process as for the “Mexican 
hat”. The energy contours however show the same kind 
of physics as a function of energy. Fig. 7(c) shows these 
contours as a function of energy for the valence band: 
again, at high negative energies, the band structure ex¬ 
hibits one unique and continuous energy contour. While 
rising the Fermi level towards the gap, the contour breaks 
up into three pockets, revealing the presence of a Lifshitz 
transition [2]. This time, however, the energy scale is 
bigger: the band gap stretches the three outer peaks and 
the difference between the three maxima at the top of 
the valence band and the minimum in the center is now 
around clt = 3 meV, as illustrated in Fig. 7(b). One 
can therefore conclude from this estimate that the ex¬ 
perimental observation of the Lifshitz transition should 
be facilitated by the layer-asymmetry induced band gap 
and that, the wider the gap, the more enhanced the vis¬ 
ibility of this effect. In Fig. 6 , we additionally show the 
density of states of such a system. This is highlighted 
with the red solid curve, which is to be compared with 
the black dashed curve corresponding to the case where 
no trigonal warping was included {vs = 0 ). 

This new topology influences the quantum Hall effect 
as well. Instead of the extra fourfold symmetry arising 
from the four Dirac cones, the triplet at the top of the 
valence band gives rise to a threefold degeneracy [ 21 ]. 
Because of the induced band gap, the valley degeneracy 
is lifted and therefore the plateaus that would survive 
down to the lowest magnetic fields would be v = ±6 and, 
because of exchange interaction, v = ±3. The full LL 
spectrum will be shown and discussed in more details in 
Section IIIB. 
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FIG. 7. (a) Band structure obtained while applying an interlayer asymmetry u = 100 meV on a bilayer graphene and taking 
into account 73 : close to the gap, the trigonal distortion is still visible, (b) Zooming at the top of the valence band, one can 
clearly observe that there are three outer maxima and a minimum in the center. The energy between them is around 3 meV for 
u — 100 meV. (c) Taking constant energy cuts in (b), one can observe that, again, as a function of energy, the Fermi contour 
gets broken, this time into three pockets: this illustrates the Lifshitz transition in a gapped bilayer graphene system. 


In this section, we pointed out how the low energy dis¬ 
persion of bilayer graphene can be manipulated by ex¬ 
ternal parameters. Effects on the Lifshitz transition are 
strong and can in principle be experimentally observed by 
studying the degeneracies of the Landau level spectrum. 
In the next Section, we will illustrate how the proposed 
effects are seen in experiments by focusing on the effect 
of the interlayer asymmetry. 

III. PROBING THE LIFSHITZ TRANSITION IN 
GAPPED BILAYER GRAPHENE 

As mentioned in the previous part, a way to tune and 
observe the Lifshitz transition in experiments is to apply 
an interlayer asymmetry to the bilayer graphene flake. 
To do so, one needs to apply a vertical displacement field 
D using a combination of top- and backgate. Encapsu¬ 
lating the bilayer graphene is therefore required. As the 
results which will be presented in this section have been 
obtained on a device fabricated from the dry transfer 
technique [38], we will first shortly describe this fabrica¬ 
tion process, but the reader should keep in mind that re¬ 
cent improvements may lead to even better performances 
[39, 40]. The electrostatics of the device will be explained 
as well as the gap properties. In a second part, we will 
reveal how, using the quantum Hall effect as a tool, one 
can identify the presence of the Lifshitz transition via the 
study of the different LL degeneracies. 

A. Inducing a gap in bilayer graphene 

The device fabrication starts with the exfoliation of 
h-BN on Si/SiC) 2 - Once a good candidate flake is iden¬ 
tified and its cleanliness is confirmed by Atomic Force 
Microscopy (AFM), a bilayer graphene flake can be trans¬ 
ferred on top. The transfer is done, as mentioned earlier, 
using the dry transfer technique: graphene is exfoliated 


on top of a ^Si-chip covered with a double-layer resist (in 
our case PVA + PMMA [50]), the first layer being water- 
soluble while the second one is a standard polymer; once 
the bilayer nature of the flake is confirmed by Raman 
spectroscopy and its cleanliness checked with AFM, the 
chip can be put in water: the sacrificial layer is dissolved 
and we are left with the graphene on PMMA floating 
at the surface. This floating polymer piece can then be 
“fished”, heated up and brought in contact with the h- 
BN using a micromanipulator. The flake is then etched 
and contacted using standard electron beam lithography. 
However, at the end of this step, the graphene surface 
is usually covered with resist residues. Different clean¬ 
ing techniques have been established: thermal annealing 
[38, 51], current annealing [52] or mechanical cleaning 
[53] . The presented device has been cleaned using the last 
technique. Once no resist residue can be identified on the 
surface, the device is encapsulated: a second h-BN flake 
is transferred on top in the same way as the graphene, 
and a top-gate patterned with electron beam lithogra¬ 
phy in the middle. The resulting device is sketched in 
Fig. 8(a) and shown in Fig. 8(b). 

Such a device defines three regions in series: the two 
outer ones, close to the Ohmic contacts, are bilayer 
graphene regions which can only be tuned using the 
doped silicon backgate, while the central area, under the 
top-gate, can be tuned combining the top- and backgate 
voltages (the carrier types - n or p - belonging to this 
area are labeled with a prime in Fig. 9). The use of 
these two gates also allows the application of an inter¬ 
layer asymmetry u and therefore allows opening a band 
gap in this region. 

In Fig. 9, the electrical characteristics of such a de¬ 
vice are shown: the conductance is measured as a func¬ 
tion of top- and backgate voltages. Two horizontal lines 
of low conductance are seen: they represent the charge 
neutrality in the two untop-gated regions: at the top of 
these lines the two external regions of the device are n- 
doped and while below, they are p-doped. Along the 
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FIG. 8. (a) Sketch of the device. A contacted bilayer graphene 
flake is sandwiched between two h-BN flakes. The whole de¬ 
vice can be tuned using the Si backgate and the central re¬ 
gion can also be independently controlled using the central 
top-gate. Combining the actions of the top- and backgate, a 
band gap is open in the central region, (b) Optical Microscope 
image of the real device. 


diagonal, another low conductance line is seen: this line, 
which defines the displacement field axis D (e.g. the line 
along which the asymmetry between top- and bottom- 
layer is increased), illustrates the situation where the 
Fermi energy lies in the gap of the central region. At 
the left of the D-axis, this area is p-doped while at the 
right it is n-doped. Along this axis however, the band 
gap size is increased, from no gap at the crossing of 
the two displayed axes, D and n gate: to a gap size of 
u ~ 80 meV at D = —0.9 V/nm (intersection of the D- 
axis and Vbg = —50 V, highlighted with a star in Fig. 9). 

The estimation of the gap size is done using a self- 
consistent procedure, as presented in Refs. [46, 49]. This 
allows to take into account the charge redistribution be¬ 
tween the layers as a function of changing gap size, which 
leads to a different screening of the applied D-field. The 
good agreement between theoretical prediction and the 
induced gap size in the fabricated devices was confirmed 
in the past by ARPES experiments [54] . From transport 
experiments, the gap size should nevertheless be accessi¬ 
ble by performing temperature dependent measurements 
(transport through the gap should obey an Arrhenius 
law). However, as demonstrated in [45, 55-57], transport 
through the gap is actually dominated by hopping pro¬ 
cesses, which obscure the true gap size. There is therefore 
no experimental way to directly confirm the estimated 
gap size in the device. We can nevertheless confirm that 
our estimation of u = 80 meV is not too far from reality, 
as confirmed both by the quantum Hall effect measure¬ 
ments which will be shown in the next section and by 
other experiments carried out on the same device, namely 
the observation and characterization of Fabry-Perot in¬ 
terference [58], the periodicity of which depends on the 



FIG. 9. Conductance map measured as a function of the 
voltages applied to the backgate, Vbg, and the one applied 
to the top-gate, Vbg- The star shows the displacement field 
value D = —0.9 V/nm, which is commented on in the main 
text. 


gap size. 

The strength of such an encapsulated geometry lies in 
the following arrangement: the h-BN material is so ro¬ 
bust against electric fields that it is experimentally possi¬ 
ble to open such large gaps, by applying such large volt¬ 
ages. This is for example not possible with suspended 
devices, which would collapse due to electrostatic forces. 
In this case, the flake remains in a good state and we are 
able to investigate this high displacement field regime. 
Every requirement is met to probe the Lifshitz transi¬ 
tion. 


B. Probing the Lifshitz transition using the 
quantum Hall effect 

As mentioned in Section IIC 2 , the influence of 73 on 
the band structure leads to a situation amenable for the 
observation of a Lifshitz transition: from one unique and 
continuous Fermi contour at positive or negative energy 
to a broken contour with three separate pockets close to 
the edges of the bands. This leads to consequences for 
the quantization observed in the quantum Hall regime: 
at low magnetic-field a sixfold degeneracy should be ob¬ 
served (three without spin degeneracy). Thus a way to 
confirm the presence of the Lifshitz transition is to per¬ 
form quantum Hall measurements on the previously pre¬ 
sented device. The results discussed in this section can 
be found in more detail in reference [ 21 ]. 

Fig. 10(a) shows the conductance map taken at 6 T. 
This map is not qualitatively different from the map 
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FIG. 10. (a) Conductance map at 6 T measured as a func¬ 
tion of the voltages applied to the backgate, Vbg , and the 
one applied to the top-gate, Vtg- (b) Measured normalized 
transconductance map: a number of lines are revealed at the 
transition between quantum Hall plateaus. 


shown in Fig. 9. However the conductance does get quan¬ 
tized and this becomes more visible when recording the 
normalized transconductance signal dG/dV^Q\ two tri¬ 
angular regions are seen, exhibiting a series of strong lines 
running parallel to the displacement field axis. These 
lines delimit a regime in which the density in the two 
outer regions is bigger than the density in the central one, 
e.g. more edge modes are present in these two regions 
(v > z/). Hence, the measured conductance, which re¬ 
flects the amount of edge modes which can travel through 
the whole device, will be given by G = e 2 v'/h [59-62]: 
in this regime, the measured conductance, and therefore 
transconductance, reflects the quantum Hall effect in the 
dual-gated region. One can then interpret each of these 
lines as the transition between plateaus of quantized con¬ 
ductance belonging to the gapped central area. 

To make this more obvious, we investigate a conduc¬ 
tance cut at high displacement fields (Vbg = —1 61 V). 
Fig. 11(a) shows in red the corresponding conductance 
measured close to 6 T: it exhibits plateaus at multiples 
of e 2 //i, indicating that all the degeneracies are lifted at 
high magnetic fields. This reveals the presence of strong 
electron-electron interaction and therefore the good qual¬ 
ity of the sample. However, a similar behavior is not ob¬ 
served on the electron side, where the LLs remain 4-fold 
degenerate. In the following, we will therefore focus on 
the hole side of the dual-gated region. 

In Fig. 11(b), the Landau level spectrum, measured 
between 0 and 6 T at Vbg = —61 V, is shown, allow¬ 
ing one to follow the evolution of our broken-symmetry 
states as a function of magnetic field. This time, the nor¬ 
malized transconductance is displayed and the plateaus 
are therefore represented by the zero-transconductance 
regions. Going from high magnetic fields to lower values, 
one notices a first unexpected feature: following filling 
factors u' = —4 and u' = —5, it appears that, close to 
5 T, these two quantum Hall states merge. Going to the 
lowest magnetic fields (B « 2.5 T) where plateaus are 
still present, two strong plateaus are seen. Taking cuts 
around this magnetic field value (yellow dashed line) and 



0 2 4 6 
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FIG. 11. (a) Conductance cuts taken along Vbg = —61 V, 
measured sweeping Vtg ~ n ga te, at two different magnetic 
field ranges: the measurements reveal a lifting of all the de¬ 
generacies at high magnetic fields and a threefold degeneracy 
at lower fields, (b) Measured Landau level spectrum, show¬ 
ing the evolution of the red conductance cuts from (a) as a 
function of magnetic field, (c) Corresponding calculated Lan¬ 
dau level spectrum, obtained for a gapped trigonally-warped 
bilayer graphene system, with an interlayer asymmetry set to 
u = 82 meV. (d) Valence band of the bilayer graphene, close 
to the top of the valence. Some representative energy cuts are 
taken from it to highlight the possible orbital degeneracies of 
the system. 


displaying the conductance, as done in Fig. 11(a) (orange 
cuts), reveals that these two plateaus are the quantum 
Hall states z/ = —3 and z/ = —6. This threefold sym¬ 
metry indicates the presence of three spin-split orbits at 
the same energy, which could be related, as mentioned 
earlier, to the presence of the Lifshitz transition in the 
gapped bilayer graphene system. 

To get a better insight into these observations, 
Fig. 11(c) shows the calculated Landau level spectrum 
corresponding to bilayer graphene, obtained by taking 
into account the presence of the gap (calculated self- 
consistently according to the D- and 5-field values) and 
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the influence of the skew interlayer hopping 73 (see B for 
more details). This calculation assumes spin-degeneracy, 
meaning that each branch in the spectrum would count 
twice if compared to the above measurement. Start¬ 
ing from low magnetic fields, this spectrum exhibits a 
threefold orbital degeneracy corresponding to the gap at 
z/ = —6 in the experiment. As mentioned earlier, this is 
due to the presence of the triplet at the top of the va¬ 
lence band, which gives rise, in magnetic field, to three 
equivalent orbits, separately confined within the three 
pockets defined by this triplet (see Fig. ll(d)i). We ten¬ 
tatively ascribe the observed gap at z/ = —3 to lifted spin 
degeneracy due to exchange interaction. While increas¬ 
ing the magnetic field, the spin-degenerate orbits grow 
bigger and start to mix in momentum space, leading to 
a splitting of the triplet, until this degeneracy is finally 
fully lifted (already at about 4.6 T in the experiment). 
At even higher fields (B = 6 T) only one unique orbit is 
left (to be compared with Fig. ll(d)iii). In-between these 
two fields, the Lifshitz transition takes place: in the band 
structure, as shown in Fig. 11 (d)ii, two contours coexist: 
the Fermi sea has, in its center, an electron-like island. 
This gives rise to two separated counter-propagating or¬ 
bits which do not mix and hence lead to an occasional 
orbital degeneracy (point encircled in Fig. 11(c): in the 
calculated spectrum, the gap of z/ = —4 closes and corre¬ 
spondingly the quantized plateau in the experiment dis¬ 
appears (dotted arrow). The position of this crossing 
corresponds fairly well to the position of the measured 
one (see yellow arrow), indicating that the value of the 
gap in the system must be quite close to the estimated 


one. The investigation of the Landau level degeneracies 
is therefore an appropriate tool to locate the presence 
of the Lifshitz transition in bilayer graphene. Beyond 
that, it is observed that odd filling factors 1 / = — 3 and 
z/ = — 5 remain quantized at the expected values as the 
Landau levels cross. Also, z/ = — 3 vanishes at magnetic 
fields slightly below the field of the Lifshitz transition. 
This behavior of the odd filling factors may be due to 
interaction effects and remains to be investigated further 
at a theoretical level. 


IV. CONCLUSION 

Bilayer graphene is a very tunable material. Small fea¬ 
tures, such as the Lifshitz transition, can be largely in¬ 
fluenced by external parameters such as strain, displace¬ 
ment fields or even magnetic fields. Bilayer graphene 
represents therefore a unique system in which the topol¬ 
ogy of the band structure can be externally influenced 
and chosen. 
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Appendix A: Bilayer graphene in quantizing 
magnetic field 


In the presence of an external, perpendicular magnetic 
field B = (0, 0, B) = V x A, we introduce canonical mo¬ 
mentum p = — iV-heA taking into account the magnetic 
vector potential A. In the Landau gauge, A = (0, Bx , 0), 
operators 7 r and 7 r^ in the Hamiltonian in Eq. (1) become 
lowering/raising operators for the Landau functions i/; n , 

[14] 


= i^vWTT)lP n+ l, 
As 

7 Tlpn = —*T—V / 2nV’r»-l; 

As 

= e^y 


\/2 ra n! A 


ex P( ~ Q X b) A b 


2A| 


Ab = el?’ 


(Al) 


where H n {x) is the n-th Hermite polynomial. Using 
these relations, one can notice that in the absence of 
the skew interlayer coupling 73 {vs = 0 ), the eigen¬ 
state of the Hamiltonian, Eq. (1), corresponding to 
the n-th Landau level, can be written in the form 
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{cl’lpn,C 2 tpn-2,C3’lpn-l,C4’lpn-l) T , where Cj, i = 1,2,3, 4, 
are some complex numbers (for the case of n = 0,1, it is 
enough to put 0 in place of any V’mco)- The coupling 73 
mixes Landau levels n and n— 3, so that the wave function 
component on each sublattice becomes a linear combina¬ 
tion of an infinite number of Landau functions ijj n . Our 
procedure to take 73 into account in the calculation of 
the Landau levels follows earlier work which investigated 
the influence of magnetic fields on the band structure of 
graphite [16]. We determine the dispersion by coupling 
n ma x Landau levels (for the spectra shown in this work, 
n ma x > 300), so that for, e.g., sublattice B2 in valley K, 
we consider n max —1 Landau functions (n = 0,..., n max — 
2). Arranging the Landau functions in the or- 
der (V’o 11 ,-ipi 1 , ipo 2 , V# 1 , i/7, ^o 2 , 4>i 2 , V’f 1 , • • • ) T in the 
valley K and (^ 2 , tpf ' 2 , V’® 1 ,-V^ 2 , V’f 2 , “V’o 11 , V>f : S 
— ipf 2 , ■ ■ ■ ) T in K' (where we used the subscripts to ex- 
plicitly denote the sublattice associated with a given < 0 n ), 
and calculating the respective matrix elements, we obtain 
the following matrix H num (in which we already included 
the interlayer asymmetry u and the strain-induced term 
w) which needs to be diagonalized numerically 
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x — vVheB , X3 = v^VheB. 


(A2) 


Appendix B: Self-consistent calculation of u 

In the presence of an external displacement field D ap¬ 
plied perpendicularly to the bilayer graphene, electrons 




FIG. 12. (a) Schematic of bilayer graphene with interlayer 
spacing d, in external displacement field D leading to the re¬ 
arrangement of the densities m and n 2 on the bottom and top 
layer, respectively. This rearrangement screens the external 
field, (b) Interlayer asymmetry calculated as a function of D. 


minimize their potential energy by rearranging them¬ 
selves between the graphene layers [46, 49]. As a result, 
the external field is screened, yielding an effective elec¬ 
tric field between the layers of magnitude E, as shown in 
Fig. 12(a). The interlayer asymmetry u can be related 
to this electric field, following 


cd 

u = ei — €2 = eEd = —[D + e(ni — 77,2)], 
e 0 


(Bl) 


where e\ (€ 2 ) is the on-site energy on the first (second) 
layer, e is the absolute value of electron charge, d is the 
interlayer spacing, cq is the permittivity of vacuum, and 
ni (77,2) is the charge density on the first (second) layer. 
Because the charge densities ni and ri 2 depend on the 
details of the band structure, which is, in turn, influenced 
by 77, Eq. (Bl) needs to be solved self-consistently. 

In the presence of a magnetic field, the difference (77,1 — 
77,2) can be determined using the capacity of a single (spin 
degenerate) Landau level, 2 eB/h, where h is the Planck 
constant, and the knowledge of the electronic eigenstates 
as obtained in A. As explained there, in the presence of 
the trigonal warping, the component c l n of the electron 
wave function of the n-th Landau level on sublattice i 
is a linear combination of the Landau level functions, 
c l n = £. The charge density difference between 

the layers is then 
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(m - n 2 ) = 2— ^(K 1 / + l C n,j| 2 - \ C n 2 j I 2 - l c n 2 | 2 )> 
n, j 

(B2) 

with the summation over n taking into account only 
those Landau levels that are filled for a specified filling 
factor v'. In the numerical calculation, we assume an 
initial value of iq n , which determines through Eq. (Bl) 
and (B2) the output asymmetry u out . We then search for 
Ui n such that the self-consistent value of u , fulfilling the 
condition u out = Ui n , is found. The number of the eigen¬ 
values included in the summation in Eq. (B2) needs to 
be varied with magnetic field, so that big enough range 
of the bilayer spectrum, ~ 1 eV away from the neutrality 
point, is included in the calculation. Because the calcula¬ 
tion is based on the single-particle Hamiltonian [Eq. (5)] 
and the spin degeneracy is implied, we are limited to even 


filling factors. At the same time, experimental features 
described in the main text are related to filling factors 
v f = — 1 to i/ = —6. Hence, we show in Fig. 12(b) the 
interlayer asymmetry u as a function of the displacement 
field D for filling factors v' = 0 and z/ = —6, to pro¬ 
vide estimates on the range of u to be expected in the 
experiment. Note that the Landau levels at the edge 
of the band, originating in the zero-energy state of gap¬ 
less BLG, are localised dominantly on one of the layers 
[14, 63]. As a result, their contribution to one of the 
charge densities ni/ri 2 is weakly dependent on u , and 
significant compared to the contribution of deeper Lan¬ 
dau levels, shared more equally between both layers. At 
small displacement fields D , this leads to large charge 
density that is placed on one layer only and cannot be 
redistributed through the self-consistent scheme, result¬ 
ing in unphysical results and lack of convergence. Hence, 
in Fig. 12(b) we do not extend the displacement axis to 
zero. Our calculation also neglects exchange effects, what 
is only justified for large u. 


